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Abstract 

We present a scenario where neutrino masses and Dark Matter are related due to a global 
U(1)b-l symmetry. Specifically we consider neutrino mass generation via the Zee-Babu two-loop 
mechanism, augmented by a scalar singlet whose VEV breaks the global U(1)b-l symmetry. In 
order to obtain a Dark Matter candidate we introduce two Standard Model singlet fermions. They 
form a Dirac particle and are stable because of a remnant Z2 symmetry. Hence, in this model the 
stability of Dark Matter follows from the global U(1)b-l symmetry. We discuss the Dark Matter 
phenomenology of the model, and compare it to similar models based on gauged U(1)b-l- We 
argue that in contrast to the gauged versions, the model based on the global symmetry does not 
suffer from severe constraints from Z' searches. 
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1. INTRODUCTION 



Observations of neutrino oscillations [1-4] imply that the Standard Model (SM) has to be 
extended in order to give mass to neutrinos. The scale for neutrino masses set by the larger of 
the two measured mass-squared differences is m v ~ 0.05 eV. Another hint for physics beyond 
the SM comes from various cosmological and astrophysical observations, which require the 
existence of Dark Matter (DM), presumably a new kind of particle beyond the known ones. 
From WMAP CMB measurements combined with other cosmological observables one obtains 
for the fraction of the Dark Matter density (Odm) to the total energy density of the Universe 
DM /* 2 = 0.1123 ±0.0035 while for normal baryonic matter one has only about 20% thereof: 
il b h 2 = 0.0226 ± 0.0005 [5]. Here h is the Hubble parameter in units of 100 km s" 1 Mpc -1 , 
with h 2 ~ 0.5. A natural candidate for a Dark Matter particle is a weakly interacting 
massive particle (WIMP), since a stable particle with mass at the "weak scale" (of order 
100 GeV) and a "typical" annihilation cross section acquires a relic abundance through 
thermal freeze-out in the early Universe with f2 DM in rough agreement with observations. 

At first sight the two scales indicated by neutrino mass and WIMP Dark Matter are 
vastly different, by about 12 orders of magnitude. Nevertheless, it is tempting to consider 
a common origin of them. In particular, it might be possible that neutrino masses are 
generated by physics at the TeV scale. In such connection between the neutrino 

mass mechanism and Dark Matter may exist, and — in the best of all worlds — both of 
them might be testable at LHC. In this Letter we present a common framework linking these 
two mysteries and offering characteristic signatures at LHC. 

There exist a number of publications which relate neutrino mass and Dark Matter; a 
few examples can be found in refs. [6-25]. We provide another interesting scenario where a 
global U(1)b-l plays a central role. For the neutrino mass mechanism we depart from the 
Zee-Babu model [26-28], where two SU(2) L singlet scalars are introduced, one single and 
one double charged, and neutrino masses are generated at two-loop level. Phenomeno logical 
studies of this model have been performed, e.g., in refs. [29-32]. We extend this model in 
a simple way by adding one scalar and two fermionic singlets. In the spirit of Majoron 
models (see e.g., [33, 34]) we impose the conservation of lepton number, or actually B — L 
(baryon minus lepton) number at the quantum level. The corresponding U(l) symmetry 
is spontaneously broken by the non-zero vacuum expectation value (VEV) of the scalar 
singlet [35]. This induces the tri-linear term in the scalar potential necessary for neutrino 
mass generation and provides a (Dirac) mass for the Dark Matter. This leads to a common 
source for lepton number breaking (and hence the generation of a Majorana neutrino mass) 
and Dark Matter. 

Any model for Dark Matter has to address the question of why the DM particle is stable 
on cosmological time scales. This is often achieved by introducing a Z2 symmetry which 
protects the DM particle, a famous example being R-parity in Supersymmetric theories. The 
adhoc postulation of such a Z2 symmetry may appear artificial, but our model provides an 
example where the Z2 symmetry emerges naturally as an unbroken remnant of some larger 
symmetry of the theory, namely the global U(1)b-l symmetry. This can be considered as a 
specific example of a more general class of models, where the dark sector respects a global 
U(l) (accidental or imposed), which gets spontaneously broken in such a way that a Z2 
symmetry remains unbroken and guarantees the stability of Dark Matter. 
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The plan of the Letter is as follows. We start by discussing the model in sec. 2. In 
sec. 2.1 we present the Zee-Babu model extended with the global U(1)b-l symmetry and 
in sec. 2.2 we discuss the scalar sector of the model. This includes a discussion of various 
consequences of the massless Goldstone boson, the Majoron, from the breaking of the global 
U(l) symmetry. We comment on implications for BBN, Higgs searches at colliders, and 
baryon genesis. In sec. 3 we introduce the DM fermions to the model. Sees. 3.1 and 3.2 
deal with the relic DM abundance and direct detection prospects, respectively. In sees. 3.3 
and 3.4 we briefly comment on indirect detection signals and on Majoron induced DM self- 
interactions, respectively. A general discussion follows in Sec. 4, where we also compare our 
model based on the global U(1)b-l to models using a gauged U(1)b-l symmetry. 



2. THE MODEL 

2.1. The Zee-Babu model and spontaneous breaking of B — L 

The scalar sector of the original Zee-Babu model [26-28] contains in addition to the 
Standard Model Higgs doublet two complex SU(2)l singlet scalars: a singly charged scalar 
h + and a doubly charged scalar k ++ . Then h + couples to left-handed lepton doublets L and 
k ++ couples to right-handed leptons £r such that the contribution to the Lagrangian is 

Aept = f tj LjC- 1 ia 2 L 3 h + + 9ij f Rl C~ l t Rl k ++ + h.c. , (1) 

where C is the charge conjugation matrix, i, j label flavour indices, and the Yukawa couplings 
/ and g are antisymmetric and symmetric, respectively. If both, h + and k ++ , are assigned 
lepton number —2, lepton number is conserved by these interactions and therefore the theory 
respects a global U(l) symmetry associated with B — L conservation. 1 In order to break 
lepton number and generate a Majorana mass term for neutrinos, usually a tri-linear term 
is introduced in the scalar potential, /j,k ++ h~h~ , which breaks B — L explicitly. Here we 
forbid such a term by demanding that B — L is an unbroken symmetry of the Lagrangian. 
Instead we introduce a complex scalar <p, which is singlet under the SM gauge group and 
has lepton number —2. Hence the VEV of (p breaks B — L spontaneously [35]. The scalar 
potential contains a term 

= X^k ++ h~h- + h.c. , (2) 
inducing the "//-term" once ip acquires its VEV, with 

W W 

Light neutrino masses are generated via a two-loop diagram, which yields 

{m v )ab = 16nfacm c g* cd I cd m d f bd , (4) 

where m c are charged lepton masses and I cd is a two-loop integral [36] . Data from the LEP 
and Tevatron colliders require that the masses of charged scalars m h and m k are typically 



Note that we only write down the leptonic part of the Lagrangian. In a complete theory baryon number 
would have to be considered as well and sphaleron processes would break B + L, while B — L is preserved. 
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larger than (9(100 GeV) [37]. Hence we can neglect the masses of charged leptons compared 
to them. In this case, one finds 

1 1 7r 2 ~ (m 



cd ~ 1 ~ T^^ll^^r 1 ) > ( 5 ) 



where M = max(mfc,m/j) and I(r) is a dimensionless function of order unity [31]. Using 
eq. (5), the light neutrino mass matrix becomes 

m -4^W' fD ' g ' D ' fT - (6) 

where the matrix De = diag(m e , m^, m T ) contains the charged-lepton masses. Due to the 
antisymmetry of /, we have detm^ = 0, and therefore, one of the neutrinos is massless as 
long as higher-order corrections are not considered. The neutrino phenomenology as well as 
other signatures of the model have been studied, e.g., in refs. [29-32]. 

Light neutrino masses are suppressed by the heavy scalar masses and proportional to 
the scale of lepton-number violation /i, set by the B — L breaking scale w, see eq. (3). 
Assuming fi ~ M ~ A, we see from eq. (6) that for m u ~ 0.1 eV the scale A of new physics 
has to be of order 1 TeV if we demand that / ~ g ~ 0.05. This scale A for generating 
neutrino masses is much below the scale of seesaw models due to the two loop suppression 
factor l/(167r 2 ) 2 ~ 10~ 4 . Hence, A is in the range of the LHC and may thus be probed 
soon. In particular, the double charged scalar provides a clean signature at colliders, via 
the decay into two like-sign leptons. Furthermore, the exchange of the charged scalars leads 
to enhanced signals in charged lepton flavour violation (/x — >■ e^y or /j, — > 3e) with good 
prospects to be observed soon [30, 31]. 



2.2. The scalar sector and the Majoron 
The full scalar potential of the model is 

Kcalar = Vp, + lA^f + A^V + f4 k++ k~~ + {l\h + h~ 

+ Ax^V) 2 + A 2 (0 f 0) 2 + Hk ++ k-) 2 + X A (h + h-) 2 

+ A 5 (^V)(0V) + \ % {^ V ){k ++ k-) + A 7 (y9»(/i + /0 

+ \ 8 (^<j))(k ++ k— ) + X 9 (^(f))(h + h-) + X l0 (k ++ k— )(h + h-), (7) 

where is given in eq. (2), denotes the Standard Model Higgs doublet, /ij are parameters 
of mass dimension one, and A« are dimensionless couplings. The A5 term in the potential 
induces mixing between if and 0. There are two massive neutral scalars in the theory, with 
propagating mass states denoted by Hi and H 2 . They are related to the real part rj of y? 
and to the real part ( of the neutral component of by 

H{\ _ ( cos a: sinaA AA , , 

H21 \— sin a cos ay I (J ' 



where the mixing angle a is given by 
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and v denotes the VEV of the Higgs doublet 0. The masses of H 1 and H 2 are 



m k 2 



= \ lW 2 + \ 2 v 2 ± \J{\iw 2 - \ 2 v 2 f + \ 2 w 2 v 2 . (10) 



The parameters pi j2 in eq. (7) can be eliminated with the help of the minimum condition 
for the potential. Therefore, the neutral scalar phenomenology depends only on three in- 
dependent parameters (in addition to the VEVs v,w), which can be chosen to be either 
(Ai,A 2 ,A 5 ) or alternatively, (m Hl ,m H2 ,a). 

Via Goldstone's theorem [38] a massless scalar has to appear in the spectrum due to the 
spontaneous breaking of the global U(l) symmetry. In our model the imaginary part of ip 
remains massless. We write _ 

V = ^(w + V + ip), (11) 

where i] = \/2Re((p) and p = \/2lm(ip) are real scalar fields and p is the Goldstone boson. 
Since the VEV of (p is responsible for lepton number breaking and the generation of a 
Majorana mass for the Dark Matter (see later) we follow the literature and call p a Majoron. 
Note that ip is a singlet under the SM gauge group and therefore there is no direct coupling 
of p to the Z boson like in triplet Majoron models [34] which are ruled out by the LEP 
measurement of the invisible Z decay width. The couplings of the Majoron p to the Higgs 
mass eigenstates Hi and H 2 are obtained from the Ai and A 5 terms of the potential (7): 

C p = ^— (m 2 Hl cos a Hi- m 2 H2 sin a H 2 ) p 2 . (12) 

The massless Majoron contributes to the relativistic energy density in the Universe, con- 
ventionally parameterized by the effective number of neutrino species. One thermalized 
scalar contributes An v = 4/7 neutrino species, and therefore a thermal scalar degree of 
freedom is consistent with bounds from Big Bang Nucleosynthesis (BBN), which are in the 
range of An^ < 1.63 [39] or An,, < 1.2 [40] at 95% CL, depending on assumptions about the 
primordial abundances. Furthermore, in our model the Majoron typically decouples from 
the plasma at temperatures above the QCD phase transition, where the effective number 
of relativistic degrees of freedom is ^ 60. Therefore, due to the entropy production at the 
QCD phase transition the Majoron abundance gets diluted and during BBN p contributes 
only with 

4 / 10. 75 x 4/3 



An u < - J « 0.06 (13) 

to the relativistic energy density, in good agreement with the above mentioned bounds. 

The prospects for discovering the charged scalars of the Zee-Babu model at the LHC have 
been discussed in detail in [31]. In particular the double charged scalar has a production 
cross section at LHC in the range of 50 to 0.1 fb for masses between 200 GeV and 1 TeV. Its 
decay into two pairs of like-sign leptons provides an essentially background free signature, see 
also [37] for a recent analysis. Possible implications of the lepton flavour violating reactions 
mediated by the double charged scalar in supernova physics have been discussed in [41]. 

In addition to the charged scalar signatures, the presence of the Majoron will modify 
Higgs physics. As pointed out in [42-44] an important feature of Majoron models are 
invisible Higgs decays H — > pp, see also [45] for a recent study. In the SM, the Higgs boson 
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dominantly decays into bottom pairs. In our model the corresponding branching ratios for 
the two Higgs mass eigenstates Hi and H 2 are weighted with the Higgs mixing angle a: 



r ( Hl -> bl) = ^ 2G l< m ^ U _ ^ ) sinV (14a) 

8tt V "l/Ti/ 



( tt ix\ 3V2G F mlm H2 ( 4m 2 . , . , 
T(H 2 ^bb) = ^ cos 2 a. (14b) 



The invisible decay modes into the Majoron p are obtained from the Lagrangian eq. (12), 

r(ifx pp) = ^^™ 3 Hl (£) 2 cos 2 a, (15a) 

r(^ 2 ^pp) = ^< (^) 2§in2 «' ( 15b ) 

and therefore [44] , 

= 1 f ^ V m 2 cot 2 a ( 1 - ^f] 1 - 250 2 cot 2 a , (16a) 



f ^) 2 tan 2 a f 1 - ^ V § « 250 ( ^) ' tan 2 a . (16b) 



T(H 2 ^bb) 12\m b J \wJ V m lf 2 , 

As we will see in the following, typically m,H 12 ~ w, and therefore, depending on the Higgs 
mixing a, we expect that one or both Higgses decay dominantly into invisible Majorons. 
In the literature, various collider signatures have been studied for the identification of an 
invisible decaying Higgs boson, for example associated production of ZH [46] and tiH [47], 
or production in weak vector boson fusion 



2.3. Baryon asymmetry 

In our model we do not attempt to explain the generation of the baryon number of 
the universe, and rely on some unspecified mechanism beyond our model. However, in 
models with low-scale lepton number violation often the problem occurs that any pre-existing 
baryon asymmetry is washed out by sphaleron processes if the lepton number violating 
reactions are in equilibrium before sphaleron freeze out. In our model lepton number is a 
global symmetry, which gets broken spontaneously by the VEV w of ip, somewhere between 
200 GeV and 2 TeV, depending on parameter values, see fig. I. Hence, at temperatures 
above w lepton number is conserved and a pre-existing asymmetry is not affected. Only 
at temperatures T < w lepton number violating processes occur via reactions involving the 
term \nwk ++ h~h~ . If sphaleron processes are still in equilibrium at that time, a pre-existing 
baryon asymmetry will be washed out. Sphaleron processes freeze out at the electro- weak 
phase transition at temperatures T ~ v. Since generically w ~ v both phase transitions 
(the B — L as well as the electro-weak ones) happen at similar temperatures. It remains 
a quantitative question of how much of the baryon asymmetry is erased by the detailed 
interplay of both phase transitions, which is beyond the scope of our Letter. Note that there 
is some parameter space where w is slightly less than v (compare fig. 1), and in this case 
B — L would be broken only after sphaleron freeze out. 
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3. DARK MATTER 



Having a global B — L symmetry which is motivated by neutrino masses, it is tempting to 
ask if it could play also a role in stabilizing Dark Matter. Therefore we introduce fermions 
TVj, which are singlets under the SM gauge group, but charged under U(1)b-l m such a way 
that the Yukawa interaction with the SM Higgs doublet <p, Lj<f)Ni (0 = ia 2 <f)*), is forbidden. 
Hence, our iVj cannot have lepton number +1 and are therefore not "right-handed neutrinos" 
in the conventional sense. Still we want the mass term for iVj to be generated by spontaneous 
lepton number breaking from the term 

C N = h tj if NjC~ l N 3 + h.c. . (17) 

This can be achieved by introducing two N fields, iV\ and N 2 -, and assigning them lepton 
numbers q\ and q 2 , respectively, such that q\ + q 2 = 2 with q 1 ^ q 2 ^ 1; for example q± = 1/2 
and q 2 = 3/2. Then no Yukawa term with the lepton doublets is allowed and eq. (17) leads 
to a mass matrix 

M "=U'o X ) W " h m * = X *T2' <18) 
and X x = h 12 . Hence, Aq and A^ 2 form a Dirac particle, with a pair of degenerate mass 
eigenstates with mass m x and opposite CP parity: 

Xi = ^(Aq + iV 2 ), X2 = ±^ Nl -N 2 ). (19) 

They are stable because of an accidental Z 2 symmetry which emerges as an unbroken rem- 
nant of the global U(\)b-l- 2 The interaction eq. (17) becomes diagonal in the Xi fields: 

Cn = \\ ^(xfC^xi + xlC'\ 2 ) + h.c. . (20) 

Note that in general global symmetries and non-gauge discrete symmetries (i.e., discrete 
symmetries which are not remnants of broken gauge symmetries) are expected to be violated 
by quantum gravity effects, e.g. [49], see [50] for a recent discussion. Planck scale suppressed 
higher dimensional operators are therefore expected in the low energy effective theory which 
would provide a mass to the Majoron [51] and might also induce DM decay in our model. 
If there is a dimension 6 operator suppressed by the Planck mass m p i inducing DM decay, 
the corresponding lifetime exceeds the age of the Universe. However, a Dark Matter decay 
operator of dimension 5 would be problematic and should therefore be strongly suppressed or 
absent, since an estimate for the lifetime yields r ~ m 2 pl /m x 10 7 s for m x ~ 100 GeV and 
there would be no DM left today. In the following we assume that Planck scale suppressed 



2 Note that there are more unbroken accidental symmetries in the Lagrangian. For example there is a Z% 
symmetry N± — > uNi,N2 — > uj 2 N2 with w 3 = 1. Another example is an additional U(l) symmetry with 
opposite charges for Ni and N 2 but all other fields uncharged. Those accidental symmetries emerge due 
to the B — L charge assignments of N\ t 2 and they ensure that no Majorana mass term is generated for 
them after B — L breaking. 
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operators induce DM decays only with lifetimes larger than the age of the Universe. We 
comment on the effects of a finite mass for the Majoron in sec. 4 

The only particle to which x% can couple is the scalar ip, with the coupling A x which is 
related to the DM mass via the VEV w, see eq. (18). Furthermore, the coupling of ip to the 
SM is provided via the Higgs portal proportional to A 5 (or the mixing angle a). Therefore, 
the two parameters A x and a will play an important role for DM phenomenology, as we are 
going to discuss in the following. 



3.1. Relic Dark Matter abundance 

The relic DM density QnMh 2 is determined by the thermal freeze-out of x% in the early 
Universe. It is roughly inversely proportional to the thermal average of the total annihilation 
cross section times the relative velocity: f2 DM ^ 2 — 3 x 10 _27 cm 3 s _1 /(°"' i; rei)- Note that the 
two DM particles xi an d Xi have identical couplings to </?. Hence the annihilation cross 
sections are the same and they will contribute in equal parts to the total DM density. Two 
DM particles Xi can annihilate into a quark-ant iquark pair, into SM gauge bosons, into the 
massive scalars H 12 , into the Zee-Babu scalars k ++ ,h + , as well as in the Majoron p. As 
examples we show the cross sections for annihilation via s-channel Hi 2 exchange into bb, 
W + W~ and k ++ k , as well as the t- and w-channel annihilation into pp: 



a bb V rel 



<7 WW"Vre\ 



&kk v re\ 



Xl S in 2 2a -^™L 
x 10247rm 



A 2 sin 2 2a gA f <l 
x 20 487rm 



x 



1 + 



A 



m v V- 



X v vc\ 

32tt 



[0 



[m 



Hi 



m 



H 2 , 



X 



m 



(m 



Hi, 
2 

Hi 



'is 



m 



m 



H 2 i 



Hi, 



X 



Is — m 



Hi, 



[s — m 



2m 2 w 



- 1 



H 2 > 



Am 



m H 2 ) cos a 



m Hi ) sin a\ 



[s — m 



Hi J 



(s - m 2 H2 ) 2 



v. 



rel 



^4 

x 15367rm 2 



2\3/2 
b) 



s - 4m 2 ,) 



s - 4m 2 ++ , 



(21a) 



(21b) 

(21c) 
(21d) 



We expanded s = 4m 2 /(l — f 2 c i/4) in v rc \ <C 1 to show explicitly the velocity suppression of 
the annihilation, however we keep s in the kinematical terms as well as in the denominators 
showing the resonant behavior of the s-channel cross section as a function of the DM mass 
m x . For simplicity we neglect here also the width of the resonances, which are however 
included in the numerical calculations presented below. In eq. (21a), and m& are the 
6-quark Yukawa coupling and mass, respectively. The annihilation cross section of DM 
into a ZZ pair can be obtained from eq. (21b) by replacing g — >■ g' = g/(\/2cos9w) and 
mw — > mz- The annihilation cross section into h + h~ can be obtained from eq. (21c) by 
replacing A6 — > A7 and m k ++ — > m h +. Note that the full annihilation cross section XX ~ > PP 
receives also contributions from an s-channel diagram, whereas eq. (2 Id) shows only the t- 
and ^-channel contribution for simplicity. This is the only case without the suppression by 
the Higgs masses 2 ■ 
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All annihilation cross sections are proportional to the coupling A x , however, in eq. (21c) 
A x appears together with the VEV w and has been absorbed into the DM mass via eq. (18). 
Annihilations into SM particles are controlled by the Higgs portal, and therefore proportional 
to the mixing parameter sin 2 2a. In contrast, annihilations into the Majoron and the Zee- 
Babu scalars survive even in the case of no mixing (sin 2a = 0, or A 5 = 0). Hence, these 
annihilation channels still provide a mechanism to generate the correct relic abundance at 
the electro-weak scale for sin 2a = 0, despite of Xi s being completely decoupled from the 
SM. 

We use the micrOMEGAs software package [52] to calculate the relic density H DM = 
Q Xl + Q X2 . For a given set of parameters, the left column of plots in fig. 1 shows the relative 
contribution of the various annihilation channels to the relic abundance. The middle and 
the right panels of the figure show the coupling A x and the singlet VEV w which are needed 
to obtain the correct relic abundance. Note that for a given DM mass m x , X x and w are 
simply related by eq. (18). We have chosen representative values for the two scalar masses 
m Hl — 300 GeV and m H2 = 120 GeV, but we have verified that our conclusions do not 
depend on this specific choice and hold within the full range of "reasonable" values for the 
Higgs masses. The couplings A M , Xs^qj^w have been set to one. Again we have checked 
that random values in the range from 0.1 to 1 give qualitative similar results. 

The upper two rows of plots in fig. 1 correspond to a relatively large Higgs mixing angle 
of sin a = 0.7 and there is a large coupling of DM to the SM. We find that depending 
on the DM mass various annihilation channels are important and the two resonances at 
m x ~ m Hl 2 /2 are clearly visible. In the lower two rows we use the same parameters but a 
small Higgs mixing sin a = 0.01. In this case the mass eigenstate Hi practically coincides 
with the singlet and therefore only the resonance corresponding to Hi exchange occurs. 
Furthermore, the coupling to SM particles is suppressed and the relic density is provided 
only by annihilations into the massless Majoron or, if kinematically accessible, into Hi and 
the Zee-Babu scalars A; ±=t , h ± . 

In fig. 1 we compare also two assumptions on the masses of the Zee-Babu scalars. For 
masses larger than the Higgs masses annihilations into A; ±=t and are sub-dominant, but if 
one of the charged scalars is lighter than one of the neutral scalars they can dominate DM 
annihilations, as visible in rows b) and d). Note also the modified shape of the curves close to 
the resonance if the dominant annihilation channel is into charged scalars. While the shape 
of the resonance of the cross section itself remains the same, the different shape for A x which 
gives the correct relic abundance follows from a different dependence of the annihilation cross 
section on A x and m x visible from eq. (21c), which shows an m x proportionality, compared 
to the typical \ x /m x proportionality of other channels. 

The middle and right columns of plots in fig. 1 show that away from the resonances the 
correct relic abundance is obtained for couplings in the perturbative range, 0.5 ^5 \ x ^3 1? 
and VEVs between 200 and 500 GeV. Close to the resonances the enhancement of the 
annihilation cross section has to be compensated by a large VEV (which corresponds to 
small A x ) in order to maintain the correct relic abundance. Note that if the symmetry of 
the model was gauged, there would be a Z' gauge boson with a mass set by w. Since Z' 
searches require the mass of such a new gauge boson to be above few TeV, e.g. [53-56], 
gauged versions of this model would be confined to the resonance regions, see ref. [17] for an 
example of such a model. It is an advantage of the global symmetry considered here, that 
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FIG. 1: Left: Relative contribution of annihilation channels to the relic DM abundance. Shown are 
annihilations into H\ and H2 ("H"), quarks ("q"), W and Z bosons ("W/Z"), Majorons ("p") an d Zee-Babu 
scalars ("h/k"). Middle and right columns show the DM scalar coupling A x and singlet VEV, respectively, 
which lead to the correct relic DM density. We assume mjix = 300 GeV, mjj a = 120 GeV, and all Zee-Babu 
scalar couplings are set to unity. Rows a) and b) are for a Higgs mixing angle sin a = 0.7, rows c) and d) 
for sin a = 0.01. Rows a) and c) are for Zee-Babu scalar masses = 300 GeV, m k ++ — 800 GeV, rows 
b) and d) for m h + = 120 GeV, m k+ + = 400 GeV. 



the breaking scale can be lower, and therefore this model does not suffer from the need of 
tuning the DM mass close to half of the mass of one of the Higgs mass states. 
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FIG. 2: Elastic scattering cross section a p of x off a proton p for = 300 GeV and mjj 2 = 120 GeV and 
mixing angles sin a = 0.7,0.2,0.05, according to the label. The singlet VEV w has been chosen in order to 
obtain the correct relic DM abundance. The left (right) panel corresponds to masses of the Zee-Babu scalars 
of k++ = 800 (400) GeV and h+ = 300 (120) GeV. Also shown are the CDMS [57] and XENON100 [58] 
exclusion limits. 



3.2. DM direct detection 

DM scattering on nuclei relevant for direct detection is mediated via t-channel exchange 
of the Higgs mass eigenstates Hi and H 2 . Hence, scattering is spin-independent. The elastic 
scattering cross section a v of \ off a proton p is obtained as 

_A|sin^2a /J 1\ 2 

°p- 4?r m ^{ m 2 Hi m 2 H J 9h p , W 

where m re d = m p m x / (m p + m x ) is the reduced mass of the DM-proton system and 



9h p 




r(p) 

-uA,s 



(23) 



see e.g., [52] where also values for fq are given. We observe from eq. (22) that the cross 
section is proportional to the Higgs mixing, since for a = the DM particle is decoupled 
from the SM and therefore the scattering cross section will vanish. Fig. 2 shows the DM 
cross section on a proton obtained with micrOMEGAs [52], where for given m# 12 ,a and 
m x the coupling A x (or equivalently the VEV w) has been chosen such that the correct relic 
abundance is obtained. Therefore the resonances in the annihilation cross section appear 
also in the scattering cross section. Although the Zee-Babu scalars k ++ and h + do not 
contribute directly to the scattering cross section, they affect the annihilation cross section 
relevant for the relic abundance, and therefore also DM-nucleus scattering depends indirectly 
on their masses, compare left and right panel of fig. 2. While part of the parameter space 
for large values of sin a is already excluded by present bounds [57-59], the cross section can 
always be suppressed by making a small. Since generically mixing should be sizable, one 
may expect observable signals in direct detection searches from this model, although they 
are not guaranteed. 
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3.3. DM indirect detection 

Let us comment briefly on possible signatures of our model for indirect searches for DM. 
The dominating branching fractions of the various DM annihilation channels are illustrated 
in fig. I. Generically our DM particles annihilate into a combination of gauge bosons and 
Higgs, leading to a rather typical mix of final annihilation products, with generic features 
common to many DM candidates produced thermally. In cases where DM annihilations into 
the Majoron dominate, visible annihilation products can be highly suppressed, which will 
lead to very small signals for indirect detection. 

If the charged scalars are not too heavy there might be a region in the parameter space 
where DM annihilates dominantly into k ±i: and/or h ± , which further will decay into charged 
leptons and neutrinos. In this case one might think that our model could explain anoma- 
lies in the cosmic electron/positron fluxes [60, 61], while respecting tight constraints for 
anti-protons. Note however, that a DM explanation of those excesses requires annihilation 
cross sections several orders of magnitude larger than the one needed for the thermal relic 
density. In our model no mechanism exists to enhance the cross section. Sommerfeld en- 
hancement, for example, would require typically light (GeV-like) mediator particles, which 
are not present in our model. Therefore, attempting to explain those anomalies within our 
framework would require an extension of the model. However, those anomalies may very 
well be of astrophysical origin. 

3.4. DM self-interactions 

As has been noted already in [33] the exchange of the massless Majoron can lead to a 
long-range force. Hence, it is important to check whether this leads to conflicts with bounds 
on DM self-interactions, see e.g. [62-64]. Due to the pseudo-scalar nature of the Majoron, 
the interaction is spin-dependent and the induced potential V(r) falls off like 1/r 3 , in analogy 
to the one-pion exchange potential in nuclear physics [33, 65]: 

3A 2 

v{r) = w^ [3(M)(f " 2) " ^ ' (24) 

where a±, Oi denote the spins of the two interacting particles. Using analogous estimations as 
in [62], we have checked that a potential of this type leads to less than one hard-scattering 
process during the age of the Universe for a typical DM particle in the halo of a Milky 
Way-like galaxy. Therefore, the Majoron induced self-interactions are in agreement with the 
corresponding bounds. 

4. DISCUSSION AND CONCLUSIONS 

We have discussed in this Letter a common framework for neutrino mass and Dark Matter 
where both are related to the breaking of a global U(l) symmetry associated with the 
conservation of B — L. From the perspective of neutrino masses, we studied a modification 
of the Zee-Babu model, where neutrinos obtain their masses at two-loop level after U{\)b-l 
breaking. The single and double charged scalars of this model provide interesting signatures 
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at collider experiments as well as in searches for charged lepton flavour violation [29-31]. 
It is then tempting to use the same global U(1)b-l symmetry in order to stabilize Dark 
Matter. To provide a Dark Matter candidate we introduced two Majorana fermions with 
non-standard B — L quantum numbers, such that a Yukawa term with the SM Higgs doublet 
is forbidden. They form a Dirac particle and are stable because of an accidental Z 2 symmetry. 
Thanks to the non-standard quantum numbers mentioned above in this framework the Z 2 
symmetry emerges as an unbroken remnant of U(l) b-l and has not to be imposed by hand. 

Due to the spontaneous breaking of the global U (1) symmetry we obtain a massless Gold- 
stone boson, the Majoron. Quantum gravity effects are expected to break global symmetries 
and hence Planck suppressed operators may provide a mass to the Majoron. As discussed 
in [51] there are certain constraints on operators from the requirement that the Majoron does 
not over-close the Universe. On the other hand, if Majoron masses are induced in the keV 
range, the Majoron could even provide a warm component to the Dark Matter [66, 67] with 
interesting signatures due to its decay [68, 69]. Here we have neglected this contribution to 
Dark Matter, and furthermore, we have assumed that possible Planck scale suppressed op- 
erators leading to DM decay maintain a lifetime of the Dark Matter particle large compared 
to the age of the Universe. 

There exists some literature on models for Dark Matter in the framework of gauged 
U(1)b-l, see e.g. [17, 18, 23, 70, 71] and it is interesting to compare with the scenario 
with a global U(1)b-l discussed here. In the gauged version the requirement of anomaly 
cancellation provides a natural motivation to introduce right-handed neutrinos while we 
had to postulate new fermions iVj. In general the right-handed neutrinos will have the usual 
Yukawa interactions with the SM Higgs doublet, and hence, an additional Z 2 symmetry has 
to be postulated if one (or more) of them should be the Dark Matter, while the protective 
symmetry arises more naturally in our case. Furthermore, a Z' gauge boson appears in the 
gauged version, which leads to additional constraints due to Z 1 searches at LEP, Tevatron, 
and LHC, e.g. [53-56]. For instance, from LEP-II ref. [53] finds w = M z ,/{2g B - L ) Z 3 TeV. 
Comparing this limit with the value of w needed to obtain the correct Dark Matter relic 
abundance in fig. 1, we find that the collider bound can be satisfied only close to the s-channel 
resonances, in agreement with [17]. We have checked that also additional contributions to 
the annihilation cross section due to Z' exchange as well as Z' — B kinetic mixing cannot 
change this result. Since the Dark Matter mass and the scalar masses are unrelated in these 
models the requirement for the resonance, ran rs 2m x , appears as an unnatural coincidence. 
A recent example for a model with neutrino masses at one loop and a gauged U(l) broken 
to a remnant Z 2 has been given in [24]. 

Before concluding let us stress that the model discussed here may serve as a specific 
example for a more general class of models, where Dark Matter is charged under some 
global U(l) symmetry and coupled to the Standard Model via Higgs-mixing of a "dark 
scalar" whose VEV breaks the symmetry. Many aspects of the phenomenology discussed 
here will apply to such models. Furthermore, for models based on B — L, the neutrino 
mass mechanism is not specific to the Zee-Babu model and one may think of various other 
possibilities to generate neutrino masses due to the U(1)b-l breaking linked to Dark Matter 
in the way described here. 
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